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The magnetization process of the spin-S' Heisenberg antiferromagnet on the kagome lattice is 
studied by the numerical-diagonalization method. Our numerical-diagonalization data for small 
finite-size clusters with S = 1, 3/2, 2, and 5/2 suggest that a magnetization plateau appears 
at one-third of the height of the saturation in the magnetization process irrespective of S. We 
discuss the S dependences of the edge fields and the width of the plateau in comparison with 
recent results obtained by real-space perturbation theory. 


Frustrated spin systems have attracted much attention 
from many condensed-matter physicists. One of the 
fascinating systems among them is the kagome-lattice 
Heisenberg antiferromagnet. Unfortunately, our 
understanding of this system is still far from complete 
in spite of many experimental and theoretical studies. 
In the S' = 1/2 system, in particular, discoveries of 
several realistic materials such as herbertsmithite,^’^ 
volborthite,^’^ and vesignieite^’® have accelerated 
theoretical studies.However, there remain some 
unresolved issues; one of them is the spin-gap problem 
of whether the spin excitation above the singlet ground 
state is gapped or gapless. 

On the other hand, fewer studies on S' > 1/2 cases 
have been carried out. As candidate S = 1 kagome-lattice 
systems, m-MPYNM-BF 4 ,®°’®i NaV 3 ( 0 H) 6 (S 04 ) 2 ,®^ 
[C 6 N 2 H 8 ][NH 4 ] 2 [Ni 3 F 6 ( 804 ) 2 ],®® and KV 3 Ge 209®4 are 
known. Theoretical studies®’®®“®® for the S = 1 case 
are also limited. Studies on the S > 1 cases have 
only started recently; candidate kagome-lattice systems 
of Cs 2 Mn 3 LiFi 2 ''^® for S = 2 and NaBa 2 Mn 3 Fii^^ 
for S = 5/2 have been reported, together with 

theoretical studies^^^'’’^ as well as an analysis based on 
the semiclassical limit.®® 

Under these circumstances, then, we are faced with 
a question: do any systematic behaviors exist in the 
spin-S' kagome-lattice Heisenberg antiferromagnet under 
magnetic fields? The purpose of this letter is to extract 
such systematic behavior of the magnetization processes 
of this model for various S by numerical-diagonalization 
calculations that are unbiased against approximations. 
With the same motivation, Zhitomirsky recently 
investigated the frustrated Heisenberg model under 
magnetic fields by real-space perturbation theory 
taking into account fluctuations around a classical 
configuration.^^ He found that in the magnetization 
process of the kagome-lattice antiferromagnet, the 
so-called uud state is stable at one-third of the height 
of the saturation, at which a magnetization plateau 
appears irrespective of the value of S. He also derived 
an expression for the I/S' expansion for both the 
edge fields of this height. The comparison between the 
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present numerical-diagonalization results and the results 
from real-space perturbation theory should contribute 
to our understanding of the frustration effect in the 
kagome-lattice antiferromagnet. 

The Hamiltonian that we study in this research is given 
by "H = Ho + HZeeman, where 

'Ho = Y,JS,-S„ ( 1 ) 

U,j) 

and 

Hzeeman = -H^SJ. (2) 

3 

Here, Si denotes the spin operator at site i, where the 
sites are the vertices of the kagome lattice. The spin 
operator satisfies Sf = S{S + 1). The sum of Ho runs 
over all the nearest-neighbor pairs in the kagome lattice. 
Energies are measured in units of J; hereafter, we set 
J = 1. The number of spin sites is denoted by W, where 
Ns/3 is an integer. We impose the periodic boundary 
condition for clusters with site Ng. 

We calculate the lowest energy of Ho in the subspace 
belonging to J2j Sj = M hy numerical diagonalizations 
based on the Lanczos algorithm and/or the Householder 
algorithm. The energy is denoted by E{Ns,M), where 
M takes an integer or a half odd integer value up 
to the saturation value Mg (= SNg). We often use 
the normalized magnetization m = M/Mg. Part of 
the Lanczos diagonalizations were carried out using the 
MPI-parallelized code, which was originally developed 
in the study of Haldane gaps.^® The usefulness of 
our program was previously confirmed in large-scale 
parallelized calculations.®®’ 

The magnetization process for a finite-size system is 
obtained by considering the magnetization increase from 
M to M -I- 1 in the field 

H = E{Ng,M + l)-E{Ng,M), (3) 

under the condition that the lowest-energy state with 
magnetization M and that with magnetization M -\- 1 
become the ground state in specific magnetic fields. 

First, let us present our results of the magnetization 
processes for 5=1, 3/2, 2, and 5/2; results are shown 
in Fig. 1. The maximum sizes of the clusters treated in 
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Fig. 1. (Color) Magnetization processes for 5 = 1, 3/2, 2, and 5/2 in (a), (b), (c), and (d), in which the maximum sizes are Ns = 24, 18, 
15, and 12, respectively. Black circles, red triangles, blue squares, green pentagons, light-blue inverted triangles, and violet diamonds 
linked by lines of the same color denote the cases of Ns = 9, 12, 15, 18, 21, and 24, respectively. 


this study are Ng = 24, 18, 15, and 12 for S' = 1, 3/2, 2, 
and 5/2, respectively. The cluster shapes we calculated 
are the same as those in Ref. 9. Note here that, the 
shapes for Ng = 9, 12, and 21 are rhombic and that 
the others are nonrhombic. In Fig. 1(a), the entire range 
of the cases up to Ns = 18 and part of the range in 
the case of Vs = 21 were already reported in Ref. 9 
and the rest of the range in the case of Vg = 21 and 
the results for Ng = 24 are additionally presented in 
the present study. For the cases of S > 1, there are no 
reports on numerical-diagonalization calculations of the 
magnetization processes to the best of our knowledge. 
The most noteworthy behavior is observed at one-third 
of the height of the saturation, where behavior similar 
to a magnetization plateau appears irrespective of the 
value of S. A detailed discussion concerning the edges 
and the width of this height will be given later. The next 
characteristic behavior is a jump near the saturation. 
Note here that all the states within the jump are 
numerically degenerate at the saturation field Hs{= 
6JS). A similar jump is known to occur in several 
cases.Regarding the existence or the absence 
of the degeneracy, this behavior is different from the 
magnetization jump observed in square-kagome-lattice 
and Cairo-pentagon-lattice antiferromagnets, and so 
on.50-53 Although this behavior occurs irrespective of the 
value of S, the skip of m at the jump gradually decreases 
as S is increased. Around m = (95' — 2)/(95) near the 


jump for 5=1, 3/2, 2, and 5/2, a region where the 
gradient of the magnetization process seems small may 
exist, although the behavior is very faint. To clarify the 
existence of the magnetization plateau at this height, 
future studies of larger systems are required. Note here 
that this height does not correspond to the m = 7/9 
plateau near the jump in the 5 = 1/2 case in Ref. 49 but 
that it corresponds to m = 5/9 for 5 = 1/2, at which 
the existence of another plateau has been pointed out for 
the 5=1/2 kagome-lattice antiferromagnet.^^’^"^ There 
appears an overhanging behavior at m = 17/21 only for 
Vs = 21 in the 5 = 1 case; this may be an artifact owing 
to a finite-size effect beacuse it is very small and beacuse 
it does not appear for Ng = 24. 

Hereafter, we focus our attention on the behavior at 
m = 1/3. 

Next, let us examine the system-size dependence of 
the edges at m = 1/3, i.e., the lower-held edge Hd and 
the higher-held edge Hc 2 ; results are shown in Fig. 2. 
In all cases of 5, the size dependences of Hd and Hc 2 
are not large. Recall here that in the case of 5 = 1/2, 
the discontinuous size dependence between Ng = 18 and 
21 is known to occur in Hc 2 f the decrease is about 
20%. In the case of 5 = 1, on the other hand, the 
decrease between Ng = 18 and 21 is about 6%, which 
is much smaller than that in the case of 5 = 1/2. It is 
therefore reasonable to use the values of Hd and Hc 2 
of the largest cluster for each 5 as substitutes for those 
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Fig. 2. System-size dependence of (a) Hc 2 and (b) Hd- Circles, 
triangles, squares, and diamonds denote the cases of S = 1, 3/2, 
2, and 5/2, respectively. 


of the infinite-size system without extrapolation when we 
carry out a detailed analysis. Note also that the small size 
dependences of Hd and Hc 2 are related to the properties 
of the m = lj3 states. After the m = lj3 states were 
studied in an analysis based on the semiclassical limit, 
it was pointed out that the m = 1/3 states reveal a 
nine-site structure in the unit cells of the spin states 
from the analysis of an effective Hamiltonian obtained 
by perturbation theory from the Ising limit,although 
Refs. 10 and 42 did not clarify the existence of the 
plateau or present estimates of Hd and H^- If the states 
form the nine-site structure, the states are energetically 
stable when iVs/9 is an integer; on the other hand, the 
energies are higher when iVs/9 is not an integer than 
when Ns/9 is an integer. This arguement suggests that 
for finite-size Hd and Hc 2 , Hd {Hc 2 ) becomes lower 
(higher) only when Ns/9 is an integer. However, such 
behavior is not observed in Fig. 2. Thus, the present 
results do not support the nine-site structure in the 
m = 1/3 states. The same situation was pointed out in 
Ref. 27 for the S =1/2 kagome-lattice antiferromagnet. 
For future studies, numerical data for Hd and Hd are 
shown in Table I together with the singlet ground-state 
energy E{Ns,0). 

Next, let us examine the 5 dependence of Hd and 
Hc 2 - We plotted Hd and Hd for the largest cluster as 
a function of l/S' the result is shown in Fig. 3 together 
with Hd and Hd for the Ns = 42 cluster in the S' = 1/2 


Table 1. Edge fields for the m = 1/3 height in the magnetization 
process of the spin-S' kagome-lattice Heisenberg antiferromagnet 
for the largest cluster treated in the present study. We also 
present the energy per site of the singlet ground state. 


s 

Ns 

Hd/J 

Hd/J 

E(Ns,0)/(NsJS'-‘) 

1 

24 

1.7502 

2.7430 

-1.4266894 

3/2 

18 

2.7016 

3.9827 

-1.2895265 

2 

15 

3.6697 

5.0085 

-1.2259126 

5/2 

12 

4.5944 

6.0308 

-1.1835511 



Fig. 3. 5-dependence of Hd and Hc 2 of the largest-size clusters. 
Circles and squares represent the results of Hd and Hc 2 for 
the largest-size clusters, respectively. Equations (4) and (5) are 
drawn as dashed and broken curves, respectively. 


case reported in Ref. 27 as the largest-cluster result for 
S = 1/2. We also draw the curves of the expressions 

Hd 
JS 
Hd 
JS 

derived in Ref. 44. One observes that the dependence 
of Hd changes between S' = 1 and S = 3/2 and that 
the numerically obtained Hd for S > 3/2 approaches 
a value of 2 with increasing S, which is the value for 
an infinite S, namely, the classical case. The agreement 
of lim 5 _,.oo Hd with the classical value suggests that 
the substitution explained in the above is reasonable. 
Note here that the 1/S dependence of the numerically 
obtained Hd is convex upward. This convex dependence 
is clearly different from Eq. (5). On the other hand, 
Hd shows an almost linear dependence on 1/S; the 
line seems to approach a value that is slightly smaller 
than 2, which is also the value for the classical case. 
The reason for the difference of this value from the 
classical value is unclear at present. The dependence of 
Hd may change above S = 5/2 if we assume a continuous 
dependence toward the classical limit. Thus, the 1/S 
dependence must be concave upward when S is large. 
This concave dependence is clearly different from Eq. (4). 
Note additionally that the almost smooth S-dependences 
of our results in Fig. 3 do not support the existence of 
the nine-site structure in the m = 1/3 state pointed out 
in Ref. 42. 

Finally, let us analyze the width of the height at 
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Fig. 4. Width A for the largest-size clusters plotted as a function 
of 1/5. The dashed line denotes Eq. (6). 


m = 1/3 in the magnetization process, where the width 
is given by A = iJc 2 — Hci\ the result is shown in Fig. 4, 
in which we also draw the line 


A _ 1 1 

7 “ 2 


( 6 ) 


obtained from Eqs. (4) and (5). The most striking feature 
is that the numerically obtained A increases as S is 
increased in spite of the fact that A is regarded as 
an energy gap in the magnetic field. Regarding the S 
dependence of the energy gap, it is well known that 
the Haldane gap of the integer-S' Heisenberg chain 
shows exponential decay with respect to The 

^-dependence of Haldane’s expression can be compared 
with reliable numerical estimates for spin-S' Haldane 
gaps.®^ The present increase in A is contradictory to 
the dependence of the Haldane gap. This feature is 
also different from Eq. (6). Note here that the present 
observation of A in units of J approaching a nonzero 
limit does not contradict the simple expectation that 
the plateau width will vanish in the classical limit 
because the plateau should be measured in comparison 
with F/g, which is linear in S. The nonzero limit of 
A in units of J is different from 1/2 in Eq. (6). One 
cannot, unfortunately, deny the possibility that A is 
overestimated in the present numerical-diagonalization 
study owing to the finite-size effect. This is a possible 
reason for the disagreement. Thus, studies tackling 
calculations of larger clusters should be carried out in 
future. Another possible reason for the disagreement is 
that the perturbation treatment around the classical 
configuration may be too rough to properly capture 
the essential quantum effect in the m = 1/3 state 
of the kagome-lattice antiferromagnet. The roughness 
may be reduced by additionally taking into account the 
effect from the spin waves; an examination along such a 
direction should be carried out in a future study. 

In summary, we have investigated the magnetization 
process of the general spin-S' Heisenberg antiferromagnet 
on the kagome lattice by the numerical-diagonalization 
method. We have found that a magnetization plateau 
appears at one-third of the height of the saturation 
even in the cases of large S. Our analysis of the 
edge fields and the width of the plateau suggests 


that the numerical-diagonalization results disagree with 
the equations obtained by real-space perturbation 
theory. The present study, based on an unbiased and 
non-perturbative theoretical method, presents significant 
information concerning general spin-S antiferromagnets 
with frustrations, which will contribute to future studies. 
In the case of S' = 1/2, the m = 1/3 state of several 
frustrated systems shows an interesting phase transition 
between the ferrimagnetic state and another state 
accompanied by a novel spin-flop phenomenon.50^53.62 
One of the cases is when the kagome lattice is distorted 
to the X type.^^ Euture studies taking this 
distortion into account could help us to more clearly 
estabish the relationship between the nine-site structure 
pointed out in Ref. 42 and unbiased numerical data. 
When the kagome lattice is spatially anisotropic, a 
non-Lieb-Mattis-type ferrimagnetic ground state of the 
5 = 1/2 system exists near the isotropic point.It 
should be examined in future what happens to these 
nontrivial phenomena when 5 > I. Further study on 
general spin-5 systems would greatly contribute to our 
understanding of the frustration effect in quantum spin 
systems. 
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